Global Phase Diagram of a One-Dimensional Driven Lattice Gas 
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We investigate the non-equilibrium stationary state of a translationally invariant one-dimensional 
driven lattice gas with short-range interactions. The phase diagram is found to exhibit a line of 
continuous transitions from a disordered phase to a phase with spontaneous symmetry breaking. At 
the phase transition the correlation length is infinite and density correlations decay algebraically. 
Depending on the parameters which define the dynamics, the transition either belongs to the univer- 
sality class of directed percolation or to a universality class of a growth model which preserves the 
local minimal height. Consequences of some mappings to other models, including a parity-conserving 
branching-annihilation process are briefly discussed. 
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The interplay of external driving fields and internal 
repulsive forces between particles can lead to interesting 
and unexpected phase transitions in the steady states 
of one-dimensional driven diffusive systems even if the 
interactions are only short-ranged jD. Generically, the 
presence of boundaries or single defects in driven sys- 
tems leads to shock waves and mutual blocking mecha- 
nisms which result in a breakdown of homogeneous par- 
ticle flow. Thus localized static inhomogeneities are re- 
sponsible for a variety of phenomena including first- and 
second-order phase transitions Q or spontaneous sym- 
metry breaking ||. These observations are of practical 
importance for the qualitative understanding of many- 
body systems in which the dynamic degrees of freedom 
reduce to effectively one dimension as e.g. in traffic flow 
kinetics of protein synthesis ||, gel-electrophoresis 
[3|, or interface growth of thin films |Q. 

Whether continuous phase transitions can occur also 
in spatially homogeneous non-equilibrium systems in one 
dimension is less well- understood Q| . In particular, there 
is a long-standing conjecture [pj that in systems with lo- 
cal interactions the steady states have rapidly decaying 
correlations and, like in 1-d equilibrium models, no phase 
transition accompanied by algebraically decaying corre- 
lations takes place. On the other hand, recent studies of 
more complicated driven systems of three or more species 
of particles in Id have demonstrated that phase separa- 
tion may take place in these models [jl0|| , thus proving the 
possibility of long-range order, but leaving open the is- 
sue of continuous phase transitions with algebraic decay 
of correlations. In the absence of a general framework 
for studying non-equilibrium phase transitions, analyz- 
ing specific models could provide useful insight in these 
complex phenomena. 

In this context, several translationally invariant one- 
dimensional growth models with local interactions which 
exhibit roughening transitions have recently been intro- 
duced. A common feature of these models is that one 
of the local transition rates which govern their dynam- 
ics is set to zero. The resulting roughening transition in 



one class of models belongs to the universality class of di- 
rected percolation [[ll| . In another class of growth models 
which preserve the local minimal height, the transition is 
found to belong to a different universality class JlJ,[l3| ■ It 
would be of great interest to put these classes of models 
within a unifying framework, so that the various types of 
transitions, the associated crossover phenomena and the 
global phase diagram could be studied. 

In this Letter we introduce a simple homogeneous 
driven Id lattice gas model with local dynamics. It ex- 
hibits a phase transition where correlations decay al- 
gebraically and which is accompanied by spontaneous 
symmetry breaking. The model can be mapped onto a 
growth model where the transition becomes a roughening 
transition. By varying the parameters which define its 
dynamics, some types of the transitions discussed above 
can be realized. The various transitions and the global 
phase diagram are studied. 

We consider a lattice gas which is an asymmetric ex- 
clusion process with next-nearest-neighbour interaction. 
Each lattice-site i S {1, 2, . . . , L} of a periodic chain may 
be either empty (0) or occupied by one particle of a sin- 
gle species, labeled A. The model evolves by random 
sequential updating. Particles hop to the right with con- 
stant attempt rate r (q) if the right nearest neighbour 
site is vacant and the nearest neighbour site at the left 
is occupied (empty). Unlike in the KLS- models 14 1, the 
left-hopping mechanism is different: A particle hops to 
the left with rate p = 1 — q — r only if the next-nearest- 
neighbour site is empty as well. The model is therefore 
defined by the transitions 



A A -> A A with rate r , 
A -> A with rate q , 
A -> A with rate p . 



(1) 



By identifying vacancies with up-spins and particles with 
down-spins, these dynamics may be interpreted as a non- 
equilibrium spin-relaxation process. The choice p = is 
a special case of the kinetic Ising models of Ref . H] , with 



1 



r = q = 1/2 corresponding to the totally asymmetric ex- 
clusion process (TASEP) §. In yet another mapping one 
obtains a growth model for a one-dimensional interface 
(see below). 

Our interest is in the stationary behavior of the half- 
filled system, i.e. the asymptotic state of the system 
reached at very large times. A thorough survey of the 
phase diagram yields as main features a phase (I) with 
spontaneously broken ^-symmetry between two antifer- 
romagnetic stationary states and a disordered phase (II) 
(Fig. As we shall argue below, the transition line sep- 
arating the two phases belongs to the universality class 
of directed percolation except for r — 0, where the uni- 
versality class is different. 
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FIG. 1. Schematic phase diagram with second order phase 
transition line (dotted curve) between the disordered phase 
77 and phase / with spontaneously broken ^-symmetry and 
non-vanishing order parameter A (cf. Eq. (^)). At r = 0, 
there is a transition at q = 1/2 to phase separation with two 
ferromagnetically ordered areas and spontaneous breaking of 
translational invariance. 

More specifically, we found that the stationary state 
can be calculated exactly along the four lines r = 0, 
q = 0, q = 1/2 and p = (Fig. |). (i) For q = 1/2, the 
system is disordered, and the stationary states are un- 
corrected product measures, (ii) For p — 0, the station- 
ary distribution is that of a one-dimensional Ising model 
. Correlations are short-ranged with divergent corre- 
lation lengths only at the extremal points q = 1 (phase 
separation into regimes with complete ferromagnetic or- 
der but opposite magnetization) and r = 1 (complete 
antifcrromagnetic order), respectively, (iii) The station- 
ary state along q = is also antifcrromagnetically or- 
dered, but at r = p = 1/2 there is an interesting phase 
transition in the dynamics of the system. Evidently, for 
small q, transitions between the two antiferromagnetic 
states A® A® A® . . . and ®A$A$A ... are possible with fi- 
nite probability, if the system is finite. However, for 
r > 1/2 the flipping time between these two states di- 
verges with a power law in system size L, whereas for 
r < 1/2 this flipping time diverges exponentially in sys- 
tem size. This is a signature for spontaneous symme- 
try breaking (and associated ergodicity breaking in the 
thermodynamic limit), even away from the line q = 0. 
(iv) Along r = the minimal height of the correspond- 
ing growth model is conserved |12(| . As in the related 



class of models of Ref. [|12|, the dynamics satisfies de- 
tailed balance with respect to an energy functional which 
is proportional to the area under the interface. The point 
p = q = 1/2 (corresponding to a change in the sign of 
the energy E) marks the transition from an antiferro- 
magnetic state to a state where complete phase ordering 
takes place and translational invariance is spontaneously 
broken. This transition is analogous to the wetting tran- 
sition of Ref. Q. 

This summary of exact results demonstrates the rich 
behavior that even rather simple homogeneous lattice 
gases may show and also indicates a certain degree of 
universality of these phenomena in ID non-equilibrium 
systems. Here, we want to discuss the behavior of the 
system as it crosses the phase transition line between 
the broken symmetry phase / and the disordered phase 
II. We shall focus on the line r = q with the limiting 
cases r = q = 1/2 (usual right hopping TASEP with mi- 
correlated disordered stationary state) and r = q = 
(left hopping TASEP with next-nearest-neighbour repul- 
sion and fully ordered stationary states). We performed 
Monte-Carlo simulations for half-filled periodic systems 
of size L = 2™, mostly with n = 10. Expectation values 
were averaged over 4000L rounds after a transient period 
of at least the same duration. 

We study the quantity 

Mt) = ] fdt^(-iy( ni {t))- (2) 
~' t=i 

where rij = corresponds to an empty site i and = 1 
to an occupied one. In the limit t — > oo, it corresponds 
to the non-conserved order parameter 2/L J2i(~^Y ( n i )■> 
which is the stationary difference in sublattice parti- 
cle densities (the 'staggered magnetization' in spin lan- 
guage). Because of ergodicity, the stationary value of 
the order parameter in a finite system vanishes by sym- 
metry. However, as a signature of spontaneously broken 
symmetry in the thermodynamic limit, one expects an 
initial decay to some quasi- stationary value A , before 
A eventually approaches zero for very long times (expo- 
nentially large in system size). On the other hand, in the 
disordered phase one expects an initially ordered state 
with A = 1 to rapidly disorder, i.e. one expects A to 
decay quickly to zero. 

A second quantity of interest is the stationary par- 
ticle current which, according to the definition (1) of 
the process, on the line r = q = (1 — p)/2 is given by 
j(q) = q(ni(l-n i+ i))-(l-2q){ (l-ni_i)(l-ni)n i+ i ). 
Clearly, j(0) = and j(l/2) = 1/8, up to a small finite- 
size correction of order 1 / L. The presence of spontaneous 
symmetry breaking suggests j = for all q < q c (up to 
exponentially small corrections in system size), since any 
finite current would lead to a transition between the two 
degenerate stationary states with A = ±Ao within a fi- 
nite time. 
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This intuitive picture is well-supported by our Monte 
Carlo simulations (Fig. |^). The current j vanishes in 
phase I and the order parameter Ao vanishes in phase 77. 
We find a phase transition point q c = 0.1515 ± 0.0005 for 
r — q, above which the current decays with a power law 



3 ~ (« - <lc) 1 



(3) 



with y w 1.7 ±0.1. Approaching the critical point q c from 
below, Ao decays with a power law 



A ~ (q c - q) 6 



(4) 



with 9 « 0.54 ± 0.04. To investigate whether this con- 
tinuous bulk phase transition is accompanied by spatial 
long-range order — as one would expect in an equilibrium 
system — we examine the stationary density correlation 
function C(k) = 4((njni+fe } — (n* ){rii + k )) which turns 
out to decay to a non-zero value below q c . At the critical 
point, correlations decay algebraically 



C(Jfe) ~ fc~ 7 . 
where 7 » 1.0 ± 0.1 (Fig. |) (l6 



(5) 
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FIG. 2. (a) Stationary current and (b) order parameter 
along the line r = q. Below: Correlation function C(fc) (c) 
for maximal distance k — L/2 as a function of r and q and (d) 
as a function of k for r — q = q c . The dashed line corresponds 
to a slope of —1. The measured data are connected by straight 
lines as a guide for the eye. 

We can gain further insight by considering the mapping 
to an interface model Jl5| which is described by height 
difference variables 1 — 2n,i and an additional stochas- 
tic variable h, representing the absolute height of the 
interface at some reference point. Each time a parti- 
cle hops to the right, the local height increases by two 
units (deposition), whereas hopping to the left describes 
a height decrease (evaporation) (Fig. |J). Thus the cur- 
rent gives the stationary growth velocity of the interface, 
while the density correlation function measures height- 
gradient correlations. Growth occurs at local minima 



with rate q, independently of the precise nature of the 
immediate environment. However, evaporation of parti- 
cles does not occur from a "flat" part of the interface: 
The corresponding process — » AA$$ is forbidden. 
On a coarse-grained scale this means that in a locally flat 
piece of the interface evaporation is not strong enough to 
create little craters which could then further grow. A 
similar situation (with different microscopic dynamics) 
was investigated by Alon ct al. [O], who found a phase 
transition between a smooth phase where no current flows 
and a rough, growing phase in the universality class of 
the KPZ-equation. The transition is related to directed 
percolation in 1+1 dimensions and is accompanied by 
spontaneous symmetry breaking in the height variable h. 




FIG. 3. Mapping between lattice gas dynamics and inter- 
face growth in 1 + 1 dimensions. A positive (negative) unit 
slope belongs to a vacancy (particle) . The interface flips (ver- 
tical arrows) correspond to particles hopping on the lattice 
(horizontal arrows). 

Here we find similar behavior which is most transpar- 
ent in the two limiting cases q = and q = 1/2, respec- 
tively. The limit q — ► 1/2 corresponds to the TASEP 
(growing, rough interface), which indeed describes inter- 
face growth in the KPZ universality class In the 
limit q — > 0, there is no current and one has spontaneous 
symmetry breaking between (macroscopically) flat inter- 
faces on an even or odd height level, respectively. We 
stress, however, that spontaneous symmetry breaking oc- 
curs already on the level of the particle description, i.e. 
without reference to the extra height variable. Assuming 
universality, one expects flT|| the exponent y to be given 
by the critical exponent v\\ s=s 1.73 of the DP-coherence 
time p8| and also a logarithmic divergence of the inter- 
face width w = [L^ 1 J2i(K - L- 1 V,^) 2 ] 1 / 2 . This is in 
agreement with our results in Eq. ([|) and Fig. [|. Also 
the value (Eh of the order parameter exponent 9 is con- 
sistent with the result 9 = 0.55 ± 0.05 reported in Ref. 
JnJ , thus independently confirming universality. Results 
on the correlation exponent 7 have not been reported in 
earlier work. 

The transition at r = is of a different nature. Here 
the model satisfies detailed balance and the current van- 
ishes both above and below the transition. At the phase 
transition point q c — 1/2 the lattice gas is uncorrelated. 
Using the interface representation of the model one can 
show |L3| that the interface width diverges algebraically 
with an exponent 1/3 as q approaches 1/2 from below 
(Fig. |). 
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FIG. 4. Interface width close to the phase transition points 

along (a) r — q, suggesting a logarithmic-like divergence, and 

(b) r = 0, respectively. The slope of the dashed line is —1/3. 

The understanding of 9 and of the new correlation 
exponent 7 (which have no conventional interpretation 
within the framework of directed percolation), and the 
behavior of these two quantities at the transition at r = 
have to be addressed in future work. Also the behavior 
of the system away from half-filling, where preliminary 
results suggest the disappearance of phase I, is an open 
issue. Returning to our original question we conclude 
at this point that the stationary states of homogeneous 
one-dimensional lattice gas models may exhibit contin- 
uous bulk phase transitions with an algebraic decay of 
correlations even if interactions are short-ranged. In our 
model, this transition results from dynamical constraints 
which — unlike in the KLS models — lead to a competition 
between a disordering dynamics (the right-hopping pro- 
cess) and processes forcing the system into either of two 
antiferromagnetically ordered states (the restricted left 
hopping process). For sufficiently strong ordering pro- 
cesses, the stationary current ceases to flow and sponta- 
neous symmetry breaking sets in. 

It is interesting to consider yet another mapping of 
our model, obtained by mapping particles into vacan- 
cies and vice versa on one (either even or odd) sublat- 
tice. The resulting dynamics are that of a new class of 
parity-conserving (PC) branching-annihilation processes 
%AA v± 000 and A00 # AAA with no absorbing state. In 
addition to particle-parity conservation (particle number 
modulo 2), there is a U(l) symmetry which results from 
the particle number conservation of the original hop- 
ping process. Generically, one expects parity-conserving 
branching-annihilation processes not to be in the DP uni- 
versality class, but in a distinct PC universality class p9| . 
From our results it appears that, in the presence of addi- 
tional symmetries, the picture of phase transitions in ID 
branching-annihilation processes is more complicated. 
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